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[0001] METHOD AND APPARATUS FOR DETERMINING 

SIGNAL-TO-INTERFERENCE RATIO WITH REDUCED BIAS EFFECT 

[0002] CROSS REFERENCE TO RELATED APPLICATION(S) 

[0003] Tliis application claims priority from U.S. Provisional Application No. 

60/425,367, filed on November 8, 2002, which is incorporated by reference as if fiiUy set 

forth. 

[0004] FIELD OF THE INVENTION 

[0005] The invention generally relates to a method and apparatus for finding the 
signal-to-interference ratio (SIR) in a digital communication system. IVEore 
particularly, the invention relates to SIR estimation with reduced bias contribution. 

[0006] BACKGROUND 

[0007] A SIR measurement is an important metric of quality performance for digital 
communication systems. For wireless communication systems such as Third 
Generation (3G) wireless systems, SIR measurements are used in several link 
adaptation techniques such as transmit power control and adaptive modulation and 
coding. Typically, SIR measured at a receiving device is more meaningful than at a 
transmitting device because SIR measured at a receiving device directly reflects the 
quality of communicated link signals, especially in the presence of multiple access 
interference or multipath fading channel. 

[0008] By definition, a received signal consists of a desired signal and interference. 
The interference may include other signals and thermal noise at the receiving end. 
However, the receiving device does not generally have knowledge of either signal power 
or interference power so that the receiving device needs to perform estimation of both 
signal and interference power based on received signals using a blind method. A blind 
method in SIR measurement for a given received signal refers to the signal power and 
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interference power (eventually SIR) obtained only from observation samples of the 
received signal without any training sequence or any prior knowledge of the desired 
received signal and interference in the received signal. 

[0009] There exist several approaches in performing measurement of received SIR. 
In the prior art, the signal power for a given signal is estimated by averaging the 
received signal over time, and the interference power is estimated by measuring total 
power of the received signal and then subtracting the estimated signal power from the 
total power. The SIR is then determined as the ratio between the estimated signal 
power and interference power. 

[0010] The SIR estimation for a given received signal can be performed at different 
observation points of the receiver structure, such as at the receiver antenna end, at the 
input to the data demodulator, or at the output from the data demodulator. However, 
SIR estimates measured at different locations usually have different levels of accuracy 
because the signal gain or the interference amount at one measurement location is 
likely to be different from the readings at other locations. 

[0011] The main problem in measuring the SIR of data signals is that an SIR 
estimate is likely to deviate from the corresponding true SIR value. Such inaccuracies 
in SIR estimation arise due to the following two main reasons. First, a signal and its 
interference cannot be completely separated. Second, desired signals are generally 
data-modulated, so that the SIR estimation is done in a 'T)lind" way, i.e., without prior 
knowledge of the data signal. This increases uncertainties in estimating signal power. 
[0012] In many prior art systems, SIR estimation mainly relies on a mean filter to 
calculate signal and noise power, resulting in undesirably large bias contribution. 
Grenerally, SIR estimation becomes more overestimated as SIR values are smaller, due 
primarily to a larger bias contribution. 

[0013] Typically, the k-th demodulated symbol, , as an input to a demodulator 
based SIR estimator, can be represented by: 
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yk - ^k '^^k Equation (1) 

d Q 
where denotes the k-th demodulated desired QPSK signal and denotes the total 



effective interference (including residual intra-cell interference, inter-cell interference 
and background noise effects), respectively. S refers to the signal and d is the desired 



By comparing Equation (2) to the SIR definition used in 3GPP, (i.e., 
RSCP*SF/Interference), neither RSCP nor ISCP is explicitly evaluated for the 
measurement. In other words, Equation (2) expresses the SIR measurement of a 
DPCH more explicitly than the 3GPP definition. In addition, since the SIR 
measurement is carried out on the data part of the received signal, a blind estimation 
is required due to the unknown transmit data at the receiving device. The function 
"E{ }" used herein represents an operator to estimate the statistical average (or 
expected or mean) value of a variable within the brackets "{ In the context of 
probability/statistics or communication systems, it is widely conventional to use E{x} to 
define the average (expected) value of a (random) variable x. 

[0014] While the SIR definition used in 3GPP is implicitly independent of the data 
demodulator type used in the receiving device, the SIR measurement in Equation (2) is 
implemented at the demodulator output. Thus, the SIR given in Equation (2) is likely 
to be different for different demodulator types. For example, for a given received signal 
primarily corrupted by interference, the SIR measured at a conventional matched filter 
receiver is likely to be smaller than that at an advanced receiver, such as an 
interference canceller, due to reduced interference effects. Note that the SIR at the 
demodulator output is the primary determinant of communication link performance. 



signal. The SIR is then estimated in terms of the average signal power, 




effective interference power, ^ , as: 




Equation (2) 
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However, the SIR measurement on the data portion of the received signal must deal 
with the unknown transmit data. 

[0015] Figure 1 depicts a t3^ical transmitted QPSK signal constellation where 

represents the transmitted QPSK symbol energy. For wireless systems such as 3GPP 
systems, after spreading the QPSK signal, the resulting spread signal arrives through 
a radio channel at the receiver. The received signal is then processed by the 

demodulator, which provides the demodulated symbols, J^jt for k = 1 , 2 , • • • ^ N ^^^^^ , 

where N^^^^^^ is the number of sjonbols in the data burst of the received signal. 

[0016] Taking into account the fading channel impact and demodulator gain, in the 
absence of the effective interference, the t3TDical signal constellation of soft-valued 

demodulated symbols can be observed on average, as shown in Figure 2 where 
represents the m-th demodulated signal symbol. 

[0017] In the presence of interference, the t3rpical demodulator output s3niibols can 

be represented pictorially as in Figure 3. For a given transmitted s3rmbol, , its 
output S3nnbol may fall into any point in the QPSK constellation, centering around the 
associated average demodulated symbol, E j^S/ } . In this case, the blind-based average 

power estimation on the demodulator output would be performed. When the decision 
for each demodulated symbol is made as to which symbol was sent, some decision error 
may occur most likely due to the effective interference and fading channel. For 

example, as shown in Figure 3, even though S2 was actually sent for the k-th sjnnbol, 

the interference may cause the demodulator output symbol, marked by yj^ , to become 

closer to in the l^t quadrant than the actually transmitted sjnnbol, As a result, 

an incorrect decision (i.e., a decision error) on y^, may be made. The decision error is 

the main source of error that causes the average signal power estimate, and 
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consequently the SIR estimate to be overestimated. In lower SIR range (high raw BER 
range), the average signal power (or SIR) estimate is likely to be more overestimated. 
[0018] It is therefore desirable to provide a method of performing SIR estimation 
without experiencing the disadvantages of prior art methods. 

[0019] SUMlVtARY 

[0020] The present invention provides a method and apparatus for SIR estimation 
which is more accurate than the estimation obtained by prior art methods. In an 
exemplary embodiment, the present invention implements a method, which preferably 
uses the demodulator output for the SIR estimation, noting that the primary 
determinant of conmiunication link performance appears to be the SIR at the output of 
the data demodulator. The method reduces such bias contribution in SIR estimation so 
that SIR estimates are as close to true SIR as possible. 

[0021] The inventive method preferably uses both a median filter and a mean 
filter and combines outputs from the median filter and mean filter in SIR estimation. 
In addition, advantageously, a correction term as a function of the mean and median 
values is introduced to further mitigate the bias effect. 

[0022] In accordance with a preferred embodiment of the present invention, an 
estimate of SIR of baseband signals which are received and processed by a data 
demodulator to provide demodulated signals to a SIR estimator is performed. The SIR 
estimator receives the demodulated sjnnbols fi-om the data demodulator and estimates 
the average signal power of the demodulated symbols as a function of a median based 
average power value and a mean based average power value of the demodulated 

sjrmbols for each quadrant of a quadrature phase shift kejdng (QPSK) constellation. 
The function is used to determine a minimum value m between and . The SIR 

estimator estimates the average effective interference power of the demodulated 
s3rmbols and calculates the SIR by dividing the estimated average signal power of the 
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demodulated symbols by the estimated average eflfective interference power of the 
demodulated s3nnbols. The SIR estimator reduces bias effects on SIR estimation. 



[0024] A more detailed understanding of the invention may be had from the 
following description of a preferred embodiment, to be understood in conjunction with 
the accompanying drawing wherein: 

[0025] Figure 1 shows a t5T3ical signal constellation for transmitted QPSK S3nnbols; 
[0026] Figure 2 shows a t3^ical signal constellation for averaged demodulated QPSK 
S3rmbols; 

[0027] Figure 3 is a t3^ical space representation of demodulated S3nnbols in the 

presence of interference; 

[0028] Figure 4 is an exemplary functional bloclc diagram of a SIR estimator 
operating in accordance with a preferred embodiment of the present invention; 
[0029] Figure 5 is a system block diagram including the SIR estimator of Figure 4; 
and 

[0030] Figure 6 is a flow chart including method steps implemented by the SIR 
estimator of Figure 4. 



[0023] 



BRIEF DESCRIPTION OF THE DRAWINGS 



[0031] 



ACRONYIVLS 



[0032] 



3GPP: third generation partnership project 



[0033] 



BER: blocl$: error rate 



[0034] 



BPSK: binary phase shift l^eying 

DPCH: dedicated physical channel 

ISCP: interference signal code power 

MUD: multi-user detection 

PSK: phase shift Iteying 

QAM: quadrature amplitude modulation 

QPSK: quadrature phase shift keying 



[0035] 
[0036] 
[0037] 
[0038] 
[0039] 
[0040] 
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[0041] RSCP: received signal code power 

[0042] SF: spreading factor 

[0043] SIR: signal-to-interference ratio 

[0044] UE: user equipment 



[0045] DETAILED DESCRIPTION OF THE PREFERRED EMBODI]VIENT(S) 
[0046] A preferred embodiment of the present invention as described below, provides 
a novel SIR estimation process based on a data demodulator output. The present 
invention also provides an SIR estimation apparatus. By definition, the term data 
demodulator output is understood to mean the output provided at the final stage of the 
considered data demodulator. The data demodulator processes received baseband 
signals and provides soft-valued estimates of transmitted symbols. The estimated 
symbols are further processed by other receiver functions, such as a channel decoder, in 
order for the receiving device to extract the transmitted data information. 
[0047] In the context of a 3GPP system, the demodulator may be configured as a 
multi-user detection (MUD) receiver or a single-user detection (SUD) receiver, such as 
a matched filter, a Rake receiver and an equalizer. Even though BPSK (binary phase 
shift keying) and QPSK (quadrature phase shift ke3dng) modulation schemes are 
referenced in the preferred embodiments, the present invention may be applied to 
higher order modulations like 8-PSK and 16-QAM (quadrature amplitude modulation). 
[0048] The present invention estimates average signal power such that the bias effect 
is mitigated. When QPSK symbols are transmitted, the average signal power of the k- 



th demodulated signal, JVjt , may be estimated as follows: 



Equation (3) 



where s( is the ^-thdemodvilated desired QPSK signal, Qi represents the i*^ quadrant 
in the QPSK signal constellation; Nq^ represents the number of the demodulator output 
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symbols belonging to the i*^^ quadrant region after making blind based S3niibol decisions 
respectively; and ykiQi) is the k-th output symbol, which is in the i^** quadrant. 

[0049] Equation (3) is utilized to determine the mean of the demodulator output 
s3rmbols in each quadrant of the QPSK constellation. Secondly, Equation (3) 
determines the average signal power based on the magnitude of the mean signal points 
in the individual quadrant. This two-step averaging (mean) mechanism may provide a 
good estimate for the average signal power in relatively high SIR range (equivalently 
low symbol error rate). However, as mentioned previously, as the actual SIR gets 
lower, the average signal power values become biased (overestimated) due to more 
s3nnbol decision errors, leading to overestimated SIR values (see Equation (2)) as well. 
To reduce the bias effect in the signal power estimation, another statistical parameter, 
called "median" (the middle of a distribution), is utilized as will be described in detail 
hereinafter. 

[0050] The mean and median are symmetrically distributed. Accordingly, with high 
SIR values, the mean and median of the demodulator output symbols located in each 
quadrant are almost identical since the interference experienced in the individual 
quadrant can be approximated to be normally distributed in a high SIR range. 
[0051] The median is less sensitive to extreme sample values than the mean. This 
characteristic of the median may make the median based average power closer to the 
true average power than the mean based average power especially for a highly skewed 
distribution such as a Log-normal distribution, or as the SIR gets lower and the 
distribution of the demodulator output samples in each quadrant approaches the 
skewed distribution. 

[0052] The standard deviation of the median for large samples with normal 
distribution is larger than that of the mean. The median is thus more subject to 
sampling fluctuations. Thus, when the number of random samples with normal 
distribution is large, the standard deviation of the median is generally greater than 
that of the mean. 
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[0053] Taking into account the above statistical properties of the median and mean, 
the present invention determines an average signal power estimate E{ } of the 
s3rmbols/bits as a function of the minimum value between the median value and mean 
value as follows: 

2 



min 



/ . 4 J 



1 ^ 

-■Y,\median{y^{Qi)} 



4 jj 



Equation (4) 



where is the ^-th demodulated desired QPSK signal, is the k-th demodulated 

S3nnbol, denotes the quadrants i of the QPSK constellation, median {yi^{Q-)) and 
mean{y^{Q^)) denote the median and mean values, respectively, of the S3niibols/bits in 

the i'th quadrant , and min( [median value], [mean value]) represents a minimum 
value function for determining a minimum value between the median and mean 
values. ykiQi) is a complex valued symbol, 

medianiyj^ {Q- )) = median{real{yj^ {Q. ))) + j • median{}mag{y^ {Q. ))) , and similarly 
mean{y^ {Q^ )) = mean{real{yi^ {Q. ))) + j • mean{imag{yj^ {Q. ))) . That is, the average signal 

power is equal to the magnitude squared of the minimum of the median absolute and 
mean absolute averaged over all quadrants. In Equation (4), the main reason for 
finding the minimum value between the median value and mean value of the 
demodulated symbols is described below. The main reason for using the minimum 
value between the median and mean value is to reduce the bias effect. It should be 
noted that the SIR value estimated according to the present invention can not be 
greater than the SIR derived by the mean value only, because the minimum value 
between the median and mean value is the smallest one of the median value and the 
mean value. 

[0054] The selection of the median value as the minimum mitigates the bias effect on 
estimating average signal power especially in the low SIR range. On the other hand, 
selection of the mean value as the minimum compensates for the drawback in median 
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calculation, such as being subject to the sampling fluctuations. Thus, by effectively 
combining the median and mean values of the demodulator output symbols for each 
quadrant of the QPSK constellation, the estimation performance of the average signal 
power is substantially improved. Even though the minimum value between the 
median value and the mean value is referenced in the preferred embodiments, other 
combined values from the median and mean may be used to determine the average 
signal power estimate. For example, a weighted (combined) method is as follows: 

-,2 



k1= 



a 



1 

-'Y.\median(y,(Q,)} 



+ (!-«)■ 



1 ^ 



mean 



Equation (5) 



where 0<= a <=1. 

[0055] Next, the averaged effective interference power based on Equation (4) is 
estimated. From Equations (1) and (4), the averaged effective interference power may 
be expressed by: 

2' 



Equation (6) 



where nl denotes the total effective interference, Nq^ represents the number of the 

demodulator output S3mibols belonging to the i^^ quadrant region after making blind 
based symbol decisions respectively, ykiQi) is the ^-th output symbol, which is in the 

i^^ quadrant, ^^{^^t } represents the average signal amplitude estimate and , 

for / = 1, 2, 3 and 4, respectively, represents the i-th QPSK constellation signal point 

1+y -1+7 -1-y 
denoted as follows: gi=-y^>g2= ^ ^ '^'^^"Ti"' ^^^^^ ^ ^ imaginary 

number, (e.g. j = V-T). 

[0056] SIR estimation must now be performed. From Equations (4) and (6), the 
demodulator based SIR estimate can be expressed by: 



-10- 



1-2-0434. lUS - Substitute Specification 



SIR = 









2 



rr 



mm 



^ /=i _ 1=1 



Equation (7) 



[0057] This SIR estimation has been validated via Unk-level simulations, which 
show that the performance of the SIR estimation based on Equation (7) is acceptable in 
reasonable operating SIR range. But in the low SIR range, (for instance, from 5 dB to 
0 dB or below), the biased effect due to symbol errors still appears in the SIR 
estimation so that it causes the estimated SIR to deviate from true SIR. The minimxim 
value can not completely eliminate the bias effect especially in the low SIR range, 
because there may be some symbol decision errors in determining the associated 
blindly based symbol decision. In this case, some correction on Equation (7) is required 
to meet the current standards requirement of 3GPP Working Group 4 (WG4). The 
present invention (including the correction term) exceeds this requirement. 
[0058] By a heuristic approach through IMonte-Carlo simulations, a correction term 
is introduced in the numerator (signal power term) of the above equation as a function 
of the offset between the calculated median and mean values, as follows: 



SIR = 



f 

min 

V 


" J 4 

- ^ '=1 




2 

- c 


1 

— • •< 
4 


4 1 ^Qi 




2 ~ 


> 



Equation (8) 



where C = 



1 

-Y.\median{y,{Q,)} 



1 ^ 
_^ 1=1 
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[0059] The basis for using such a correction term is that in a high SIR range, the 
median value and mean value in the correction term are mostly likely close to each 
other. Therefore, the correction term can be negligible when in the range it is not 
needed, (since without it the estimated SIR is already within the accuracy 
requirement). However, as the actual SIR gets lower, the skewed distribution due to 
symbol error effect of the demodulator output samples may drive the correction term 
increasingly, since the difference between the corresponding median and mean values 
may be gradually increased. At the same time, as the estimated signal power 
(eventually SIR) gets overestimated (biased), the correction term may help to reduce 
the biased effect in the SIR estimation. 

[0060] Although the SIR measurement method described hereinbefore was derived 
under the assumption that the demodulator output is the complex valued QPSK 
symbols sequence, for a practical MUD implementation, the MUD provides real valued 
data bits sequences with each pair of two consecutive data bits, which can map to a 
complex valued symbol such as QPSK modulation in a transmitter. 
[0061] Figure 4 shows a block diagram of a SIR estimator 400 taldng input as data 
bits rather than QPSK sjnnbols, in accordance with a preferred embodiment of the 
present invention. The SIR estimator includes an input port 405 for receiving an input 
bit sequence, a hard limiter 410, a multiplier 415, a median filter 420, a first mean 
filter 425, a minimizing process block 430, a signal power process block 435, a 
correction term process block 440, a second mean filter 445, a SIR calculating process 
block 450, summers/comparators 455, 460 and a process block 465. 
[0062] Input port 405 receives a soft-valued bit sequence. The SIR estimator 400 
processes absolute values of symbols input to the SIR estimator in the form of a bit 
sequence received via input port 405. The bit sequence is routed to the hard limiter 
410 and the multiplier 415. The hard limiter 410 provides a +1 bit to multiplier 415 if 
a sofi;-valued bit is greater or equal to zero. Otherwise the hard limiter 410 provides a 
-1 bit to the multiplier 415. Multiplier 415 then multiplies each soft-valued input bit, 
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, by the corresponding hard Hmited bit, resulting in the absolute value |^ | of the 
individual input bit. 

[0063] The absolute value |^ | indicates that a hard-bit decision on each input bit, ^ , 

has been made and if the resulting bit decision turns out to be -1, then the input bit is 
180 degree phase shifted. Otherwise, the input bit remains unchanged. Thus, the 
calculation of the average signal power and interference power is based on a blind 

based bit decision. The multiplier 415 outputs the absolute value |^ | to median filter 

420, first mean filter 425 and summer/comparator 460. Based on a running number of 
samples, the median filter 420 and mean filter 425 determine the median value and 
mean value, respectively, of the absolute valued bit sequence. The outputs from the 
median filter 420 (ma) and first mean filter 425 (me) are compared in the minimizing 
process block 430 to determine a minimum value m between the median based average 
power value ma and the mean based average power value me. The correction term 
process block 440 also receives the outputs fi-om the median filter 420 (ma) and mean 
filter 425 (me) and performs the following calculation to determine a correction term C : 

C = - m^f , Equation (9) 

J 4 1 

where = ^\median{yf^{Q^)J and rn^ = T ' X |^^^^ (^'/t , yk is the k-th symbol 

in the sequence of symbols, denotes the quadrants / of the QPSK constellation, and 
median (yf^iQ^)) and mean{yf^{Qi)) denote the median and mean values, respectively, of 

the S3nnbols in the i-th quadrant . 

[0064] The output m of minimizing process block 430 is routed to the signal power 
process block 435 and to sununer/comparator 460. The summer/comparator 455 
compares the output Ps of signal power process block 435 to the correction term C, 
where: 

P, = (mf Equation (10) 
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[0065] In order to determine the average interference power, the process block 465 
first receives the output of summer/comparator 460 and performs the function 

r^-mf which extracts the interference component out of the input bit sequence. The 

second mean filter 445 receives the output of process block 465 and outputs Pn to the 
SIR calculating process block 450. The SIR is calculated by SIR calculating process 
block 450 based on the output from summer/comparator 455 and second mean filter 
445, where: 

SIR = —Z — Equation(ll) 

[0066] Figure 5 shows a system 500 which includes a demodulator 505 and a known 
sofl;-syTnbol-to-sofl;-bit mapper 510 which inputs a soft-valued bit sequence to the SIR 
estimator 400. The SIR estimator 400 may be used also for higher order modulations 
like 8-PSK, 16-QAM, and 64-QAM, if the complex valued demodulated sjnnbols are 
converted to soft-valued bits via the soft-syrobol-to-soft-bit mapper 510. 
[0067] Figure 6 is a flow chart of a method 600 implemented by the SIR estimator 
400, The SIR estimator estimates the SIR of symbols/bits by receiving the s3rmbols/bits 
(step 605), estimating the average signal power of the symbols/bits as a function of a 
median based average power value and a mean based average power value of 

the symbols/bits for each quadrant of a QPSK constellation (step 610), estimating the 
average effective interference power of the sjnnbols/bits (step 615) and calculating the 
SIR by dividing the estimated average signal power of the symbols/bits by the 
estimated average effective interference power of the symbols/bits (step 620). The 
function of the median based average power value and the mean based average 

power value is to provide a minimum value function for determining a minimum 

value m between the median based average power value and the mean based 

average power value . The average signal power of the symbols/bits is equal to the 

magnitude squared of the minimum of the absolute value of the median based average 
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power value and the absolute value of the mean based average power value 
averaged over all of the quadrants of the QPSK constellation. 

[0068] The foregoing describes a novel SIR estimator based preferably on data 
symbols. Output of the data demodulator of interest, such as Rake output or MUD 
output, is fed to the SIR estimator. As indicated above, an advantage from using a 
data demodulator output as input is that the demodulator output best and directly 
reflects quality of received data signals. Especially when SIR measurement is used for 
link control techniques like power control, a data-demodulator based SIR measurement 
as described above is highly desirable. In addition, the proposed SIR estimator is 
capable of reducing bias effects on SIR estimation, resulting in more reliable and 
accurate SIR estimation than the prior art. All such modifications and variations are 
envisaged to be within the scope of the invention. 

[0069] While the present invention has been described in terms of the preferred 
embodiment, other variations which are within the scope of the invention as outlined 

in the claims below will be apparent to those sldlled in the art. 

* * * 
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